Optimal Lewenstein-Sanpera decomposition of two-qubit states 
using Semidefinite Programming 



On 
O 
O 
(N 

Oh' 

<D 

in 

(N 



9 L,: 



> 

On 
On 

in 

ON 

o 

On 
O 



X 



Guo Chuan Thiang, 1 ' 2 Philippe Raynal, 1 and Berthold-Georg Englert 1,2 

1 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 117543, Singapore 
Department of Physics, National University of Singapore, 2 Science Drive 3, 117542, Singapore 

(Dated: September 25, 2009) 

We use the language of semidefinite programming and duality to derive necessary and sufficient 
conditions for the optimal Lewenstein-Sanpera Decomposition (LSD) of 2-qubit states. We first 
provide a simple and natural derivation of the Wellens-Kus equations for full-rank states. Then, we 
obtain a set of necessary and sufficient conditions for the optimal decomposition of rank-3 states. 
This closes the gap between the full-rank case, where optimality conditions are given by the Wellens- 
Kus equations, and the rank-2 case, where the optimal decomposition is analytically known. We 
also give an analytic expression for the optimal LSD of a special class of rank-3 states. Finally, 
our formulation ensures efficient numerical procedures to return the optimal LSD for any arbitrary 
2-qubit state. 



I. INTRODUCTION 

Quantum entanglement is perhaps the most striking 
phenomenon associated with quantum systems. Once 
seen as "evidence" for the alleged incompleteness of quan- 
tum mechanics entanglement has now found numer- 
ous applications as a resource for quantum communica- 
tion 0, [1[ , computation [J] , and cryptography [j| . 

The characterization of entanglement has become a key 
area of quantum information theory. Various schemes to 
quantify entanglement have been proposed (for a review, 
see @). A particularly interesting approach is offered 
by the Lewenstein-Sanpera decomposition (LSD) Q of a 
composite quantum state, which comprises a convex sum 
of a separable state and an entangled state. 

Now for any 2-qubit system, there is a unique optimal 
LSD. This optimal decomposition has a separable part 
with maximal weight, and the entangled part is a pure 
state. The weight of the pure state in this decomposition 
multiplied by its concurrence 0, 0] provides a measure of 
entanglement for the 2-qubit state 0, E3] • 

Analytical expressions for the optimal LSD of some 
special cases were found in [ll|, these include the rank- 
2 states, the self-transposed states, and the generalized 
Werner states. Recently, a pair of coupled nonlinear 
equations for finding the optimal LSD of full-rank states 
was obtained by Wellens and Kus However, an an- 
alytic solution to these equations is only available in the 
case where the separable part in the optimal LSD has full 
rank. 

As noticed in [f2j , the problem of finding the optimal 
LSD can be in some cases formulated as a SemiDefinite 
Program (SDP). In the present paper, we systematically 
exploit this connection for 2-qubit states. We first red- 
erive the Wellens-Kus equations for full-rank states in 
a particularly transparent manner. The SDP formula- 
tion also enables us to efficiently compute the optimal 
decomposition by numerical means. We then extend our 
analysis to rank-3 states, and obtain necessary and suf- 
ficient optimality conditions. With the optimal LSDs of 



rank-2 states already known [13j , this completes the char- 
acterization of optimal LSDs for 2-qubit states. We also 
obtain analytically the optimal LSD for the class of rank- 
3 states that are orthogonal to a product state and have 
a separable part of rank 3. For such states, the pure state 
in the optimal LSD is maximally entangled. This is sim- 
ilar to the full-rank case where the separable part is full 
rank. There, the nonseparable pure state is maximally 
entangled too [l4| . 



II. LEWENSTEIN-SANPERA 
DECOMPOSITIONS 

The construction of LSDs hinges on the fact that the 
set of separable states is convex. Any composite system 
can be written as a convex sum of a separable state p sep 
and an entangled state p D nt- Information about nonsep- 
arability is then contained in p en t] for example, the state 
p is nonseparable if p en t does not vanish, and only then. 

A simple dimensional argument leads to the impor- 
tant consequence that for 2-qubit states, p G nt is just a 
pure state. In general, there is a continuum of LSDs, 
p = Xp sep + (1 — A)/9 purc , for a given state. Among these 
is the optimal LSD, 



p = Sg scp + (1 - S)g pu 



S = 



<{A}, (1) 



where S is the degree of separability of p. Throughout 
this paper, we will use calligraphic font to refer to quan- 
tities that are optimal. 

When p has full rank, g sep is either full-rank or rank-3. 
In the latter situation, we denote its null eigenstate by p\ . 
Let us also introduce gj ep , the partial transpose with re- 
spect to the first qubit of g SC p- Then the barely-separable 
property of g sep [LJ] says that g^L has a zero eigenvalue, 
whose corresponding null eigenstate shall be denoted by 
P2 ■ We quote the following results from the Wellens-Kus 
paper [10| . with slight modifications to their notation. 

In the optimal LSD of a full-rank state, g purc is an 
eigenstate of pp\ + pj 1 , p, > , with a nonpositive eigen- 
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value, 

3a, [I > (fipi + P2 1 )0pure = -agpure, (2) 

with fipi = if Q sep has full rank. This is accompanied 
by the eigenstate equation for p 2 , 

(p - (1 - S)g purc ) Tl p 2 = 0. (3) 

Equations ([2]) and ([3]) are the Wellens-Kus equations. 
In general, there may be several solutions to these cou- 
pled eigenvalue equations. However, consistent with the 
uniqueness of the optimal LSD, there is only one with 
p,a > that gives a positive and separable, and thus 
permissible, p scp . 

The original proofs of these assertions, as well as the 
sufficiency of these equations, involve considerable tech- 
nical detail. The aim of the present paper is to present an 
alternative derivation, and to generalize these equations 
to the reduced-rank case. 



III. SEMIDEFINITE PROGRAMMING 

In semidefmite programming [15], a linear objective 
function is minimized subject to the constraint that 
an affine combination of hermitian matrices is positive 
semidefmite. We now briefly review some important fea- 
tures of SDR 



A. The primal semidefmite program 

In its canonical form, the primal semidefmite program 
is formally stated as: 

minimize c T x 

subject to F(x) > 0, ^ ' 

where F(x) = F Q + Y%Li x i F i and ^ e R ™- The in P uts 
for the primal problem are (i) the vector c € K m char- 
acterizing the objective function, and (ii) the m + 1 her- 
mitian matrices Fq,F\,..., F m <= TL n defining the linear 
matrix inequality, where TL n is the space ofnxn Hermi- 
tian matrices. The primal problem is strictly feasible if 
there exists x such that F(x) > 0. The primal optimal 
value is p* = inf{c T x | F(x) > 0}, and we denote the 
primal optimal set by 

X opt = {x | F(x) > and c T x = p*}. (5) 

B. The dual semidefmite program 

The dual problem associated with is 
maximize — tr{ FqZ} 

subject to tr{F,Z} = Cj,i= 1, . . . , m, (6) 
Z > 0. 



The dual variable Z = G TCI is subject to m equality 
constraints, defined by the F^s and CjS specified in the 
primal program, in addition to a condition of nonnega- 
tivity. The dual problem is strictly feasible if there exists 
Z > satisfying the dual constraints. The dual optimal 
value is d* = sup { - tr{F Q Z} \Z > 0, tr{F t Z} = c t Vi}, 
while the dual optimal set is 

Z pt = {Z > | tr{F t Z} = Ci Vi, -tr{F Q Z} = d*}. (7) 

One also has the hierarchy —tr{FoZ} < d* < p* < c T x, 
meaning that the dual objective yields lower bounds 
on the optimal primal value, while the primal objective 
yields upper bounds on the optimal dual value. 

C. Complementary slackness condition 

An important quantity to consider is the duality gap 
c^-x + tr{F()Z} = tr{F(x)Z}, which is a nonnegative 
quantity linear in x and Z. The equality d* = p* holds 
(no duality gap) if either the primal or the dual problem 
is strictly feasible. If both are strictly feasible, the optimal 
sets X op t and Z opt are nonempty, and there exist feasi- 
ble pairs of x and Z with p* = c^x = — tr{FoZ} = d* , 
so that F(x)Z = 0. This is the complementary slack- 
ness condition, stating that the ranges of the nonnega- 
tive matrices F(x) and Z are orthogonal. Under strict 
primal and dual feasibility, one then has necessary and 
sufficient optimality conditions for the semidefmite pro- 
gram: a feasible x is optimal if and only if there exists a 
Z such that 

F(x) > 0, Z > 0, 

tr{FiZ} = a, i=l,...,m, (8) 
F(x)Z = 0. 

The above equations provide algebraic expressions that 
the optimal x and Z must satisfy. We will see in the 
next section that these conditions lead to the Wcllcns- 
Kus equations. 

IV. DERIVATION OF THE WELLENS-KUS 
EQUATIONS 

Let us notice that we have an optimization prob- 
lem, in which we wish to minimize a scalar function 
— A = — ti'lApsop} of some variables subject to a set of 
constraints. Firstly, we require p sep and p pU re in a LSD 
to be positive semidefmite. Next, the Peres-Horodecki 
criterion [lj| [ijj tells us that a 2-qubit state is separable 
if and only if its partial transpose is positive. The crucial 
point here is that the separability constraint has become 
a positivity constraint, ensuring that the optimal LSD 
problem for 2-qubit states can be formulated as a SDP. 
We will proceed to show this explicitly. For simplicity, 
we only consider full-rank states in this section. The case 
of reduced-rank states will be considered in the following 
section. 
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A. Optimal LSD as a semidefinite program 

1. The primal problem 

We use a and f to denote the Pauli operators in the 
first and second qubit space, respectively. It will be con- 
venient to use the magic basis, introduced by Hill and 
Wootters [8, 9], in which the Pauli operators are repre- 
sented by imaginary antisymmetric 4x4 matrices while 
their products are represented by real, symmetric matri- 
ces. Partial transposition in the first qubit is effected by 
a —> —a, t — > t. 

Our basis {Ei : i = 1, . . . , 16} for 4x4 hermitian 
operators comprises the sixteen combinations of the Pauli 
operators and the identity, <7jTj, where i, j = 0, 1, 2, 3 and 
(jo = to = fl-4- These are traceless (except E\ — I4) and 
mutually orthogonal, i.e., tr{EiEj} — iSij. A LSD of a 
state p can be written as 



Next, we establish strict primal feasibility. For this, 



P — Ap S cp + (1 — A)p purc = p S ep + f>pu 



(9) 



where the weights A and 1 — A have been absorbed into 
/5 scp = Vsop and p pure = (1 - A)p puro . In this nota- 
tion, we have the parameterization p. 



scp 



• E, where 

x 1 ' = (X,x 2 ,...,xie) GM 16 . 

In the search for the optimal LSD, we comb through 
the possible p sep s via choices of x, but these choices are 
not arbitrary. To ensure a valid decomposition in the 
first place, we must enforce three constraints, 



(i) positivity of p sep 

(ii) separability of p scp 
(hi) positivity of pp Ure 



(10) 



which we merge into a single inequality of a 12 x 12 ma- 
trix: 



Pscp 

P s Tl 







so p 

p 



> 0. 



Next, we introduce 16 block-diagonal 12 x 12 
mitian matrices Fi associated with the EiS, 
fined by F t = |diag(^, Ej 1 , -£?,-),» = 1, . . 
as well as Fo = diag(0, 0, p). In terms of the 
the inequality constraint in Eq. (jll[) can be 
pressed ' r " ' ' ' n ^ 



as F(x) =F Q + x i F * > °- Finally, 

0I6 



(11) 



her- 
de- 
,16, 
F iS , 
ex- 
let 



= (-1,0, . . . ,0) G R lb , so that c l x = -A. Maximiz- 
ing A to obtain the optimal LSD is then equivalent to 
minimizing c T x. 

With these specifications, we have rephrased the opti- 
mal LSD problem as a SDP in the form of (j4]). One can 
then efficiently compute the optimal LSD of a given 2- 
qubit state using well-established algorithms for solving 
SDPs. For instance, we have written a working routine 
using cvx version 1.2 [l8j |. which is a modeling system 
for disciplined convex programming, utilizing the open- 
source solver SDPT3 111. 



we choose 



i.e., x 



(a, 0, . . . , 0), where 



a j is a positive number smaller than the smallest eigen- 
value of p. Clearly, p sep > and p^L > 0. Furthermore, 
since p > 0, it has a spectrum with 4 positive eigen- 
values, and choosing a as described above, the differ- 
ence p — pscp = p — a^l4 is still positive definite. Thus, 
Eq. (fTT|) holds with strict inequality as required. Since 
the primal problem is strictly feasible, we conclude that 
there is no duality gap. 



2. The dual problem 

We now focus our attention on the dual problem as- 
sociated with (U]). Following ©, the dual variable Z 
is a 12 x 12 positive semidefinite matrix subject to the 
16 dual constraints tv{FiZ} = a. Since Fo and Fi are 
block-diagonal, the dual objective depends only on the 
block-diagonal entries of Z. Without loss of generality, 
we can choose Z to be block-diagonal. For convenience, 
we write 



Z = 



2i 
z 2 
z 3 



(12) 



where Z\, Z 2 and Z3 are nonnegative 4x4 matri- 
ces. With this notation, the dual objective becomes 
— tv{pZ^\. Since there is no duality gap, we have 
d* = -tr{pZ 3 } =p* = -S. 

The dual problem is strictly feasible too: choose 
Z = diag(l4, 14,314) > 0, and check that all the con- 
straints are indeed fulfilled. The first dual constraint 
tr{iqZ} = — 1 is satisfied and the 2nd to 16th dual con- 



straints tv{FiZ} 
by construction. 



hold, since Ei and E; 1 are traceless 



B. Equivalence of complementary slackness 
condition and Wellens— Kus equations 

With strict primal and dual feasibility, we now have 
necessary and sufficient optimality conditions as a con- 
sequence of the complementary slackness condition 
In the present context, conditions ([5]) translate into the 
following statement. The primal variable g sep is optimal 
if and only if there exists a Z such that 



= 0. 





Pscp 







~z x 





" 


(I) 





^ P 







z 2 










f?pure 










z 3 



(II) 



e 



scp 

Qp 



>o, 



4 



(III) 



Zi 
Z 2 
z 3 



>o, 



tr{F,Z} = a 
i = l,..., 16. 



(13) 



Here, £> S ep> £?pure, and Z refer to the optimal variables. 
Let us digest this information. (I) is a set of three eigen- 
state equations from the slackness condition that deter- 
mines the matrices Z%, Z 2 and Z 3 . (II) is the primal 
constraint and simply reiterates that we have a valid de- 
composition in the first place. (Ill) is the set of dual 
constraints, which we will utilize to express Z 3 in terms 
of Zi and Z 2 . Notice that the i^s are composed of blocks 
of EiS, the 16 orthogonal basis matrices for the space of 
4x4 hermitian matrices. In fact, the 16 dual constraints 
trji^Z} = a are really statements about the 16 compo- 
nents of the operator Z\ + Zj 1 — Z 3 in the "directions" 
of Ei. Specifically, the ith dual constraint reads 



Ei 
Eji 


1 

I 






-Ei 





Zi 





" 









z 2 





j 










z 3 





tr^Z,} 
tr{E i (Z 1 +2^ 1 -Z 3 )} = 



rtr{^Z 3 } 



(14) 



where we used the identity tr{Ej 1 Z 2 } = trj^Zj 1 }. 
Since any hermitian operator can be written as 
H =\ Y%Li E i ^{E l H}, we arrive at 



Z 3 =Zx 



Z Tl 

^2 



(15) 



We are now ready to state the Wellens-Kus equations. 
The third block equation in (II) states, using Eq. (|15p . 



rank. Then we must have Z 2 = and Eq. (| 16[) becomes 
Zigpuro = — £> pure . But Zi is assuredly nonnegative by 
(III), so £> P urc must vanish, which is to say, p was sepa- 
rable to begin with. 

Now for a nonseparable p, g^} has rank 3 so Z 2 must 
be a pure state. If in addition, g scp has full rank, Z\ must 
vanish. In this is the pure state associated 

with the negative eigenvalue of Zj 1 , which is a Bell state 
[20| . This is consistent with the observation made by 
Karnas and Lcwcnstcin in Ref. 



14] 



In passing we note that Z\ , Z 2 and Z 3 have an in- 
teresting interpretation in the language of entanglement 
witnesses. An entanglement witness W is a hermitian op- 
erator such that trjVFpsop} > for all separable states 
Pscp, but for some entangled state p en t, tr{Wp cn t} < 0. 
The dual of the optimal LSD problem for 2-qubit systems 
can be written as an optimization over a constrained set 
of entanglement witnesses [2l| , so that 



S = max{0, -min trjWp}}. 

W+ti>0 



The quantity Z\ + Zj 1 can be interpreted as the optimal 
entanglement witness W for the state p, since 



(20) 



tr{Wp sop } > V separable states p sep , 
tr{Wp} = S - 1 < 0. 



(21) 



It is optimal because tr{p sop (Zi + Zj 1 )} = 0, so 
Z\ + Zj 1 "ignores" the separable content of p, while 
maximally detecting the entangled part g pulc in accor- 
dance with Eq. (HI 



(Zi + Z 2 l )e pU re = -^pure- (16) 

This is supplemented by the second block equation in (I), 
in which we carry out the replacement g scp —> p — g pme 
to obtain 



(P - £purc) Tl Z 2 = 0. 



(17) 



Equations (fl"6"j) and ([TTJ) are the Wellens-Kus equations, 
which we restate here for easy reference: 



3a, /J, > (npi + P2 1 )Qpure = -ag pU re 

(p - (1 - S)g pulc ) Tl p 2 = 0. 



(18) 
(19) 



The first block-equation in (I) states that g scp Zi = 0, so 
Z\ is proportional to p\. Therefore, Eqs. (TT5]) and (TiTq)) 
are really the same equations, with the multiplicative fac- 
tors a and p absorbed in the normalization of Z\ and Z 2 . 
It is also clear that Eqs. (|19[) and (jTTJ) are the same equa- 
tions, with p2 and Z 2 differing only by a multiplicative 
factor. 

We remark that the barely-separable property of ^ sep 
in the optimal LSD of p can be derived as a consequence 
of this formulation. Suppose otherwise, that g^ p has full 



V. GENERALIZED WELLENS-KUS 
EQUATIONS FOR REDUCED-RANK STATES 

Since the optimal LSDs for rank-2 states are already 
known, it remains to characterize the rank-3 states to 
fully apprehend the LSD of any 2-qubit state. As a side 
result, Wellens and Kus [l0| generalized their equations 
to the reduced-rank states by treating them as the limit 
x — » of the full-rank state a^lU + (1 — x)p. However, 
their approach has the implicit assumption that Z 2 is a 
pure state, whereas it could also be of rank 2. The com- 
ponent Zj 1 in the optimal entanglement witness need 
not be the partial transpose of a pure state. As we will 
show, the SDP approach naturally takes care of this sub- 
tlety. 

Clearly, the primal problem in the previous form is 
never strictly feasible if p has rank 3. In order to uti- 
lize the complementary slackness condition, we need to 
modify the primal problem such that strict feasibility is 
restored. We denote the pure state orthogonal to p by 
7 and its concurrence by q. There will be two separate 
cases to consider: (i) 7 is entangled, and (ii) 7 is a prod- 
uct state. 
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A. 7 is an entangled state 

1. The primal problem 

We consider a parameterization in the three dimen- 
sional subspace spanned by p, which requires 3x3 = 9 
parameters. The rank-3 projector onto the orthogonal 
complement of 7 is given by P 3 — 1 4 — 7. We denote by 
I3 its restriction to its own support. In its generic form, 
7 can be written as 

7 = i(l 4 +pai -pr\ - o\T\ - qa 2 r 2 - qa 3 r 3 ), (22) 

where p = \Jl — q 2 and < q < 1. One can then con- 
struct an orthogonal basis {I\ : i — 1, ...,9} for the 
support of I3, in which r 1 = I3 and the remaining I\ 
are traceless. An explicit construction for {I\} can be 
found in [l3|]. In this basis, the parameterization for the 
(unnormalized) rank-3 state p sep becomes p sop = ^x ■ T, 
where the primal variable x — (A, x 2 , ■ ■ ■ , xg) is in R 9 . 

One can represent the I\s by 3 x 3 matrices, but their 
partial transposes Tj 1 s can be full-rank, therefore we 
need 4x4 matrices to write them. Following the same 
prescription as in the full-rank case, we express the three 
primal constraints in block diagonal form, 



"0 0" 


1 9 


Ti 





" 








If 1 





p 


6 i=l 











where the first and third blocks are 3x3 and the 
second block is 4 x 4. Analogously to the full-rank 
case, we define Fj = |diag(r i ,r^ 1 , — Ti),i — 1,...,9, 
and F = diag(0, 0,p), so that Eq. (f2"3"|) turns into 
F(x) = Fq + J2^=i x iFi — 0. Finally, we also define 
c T = (-1,0,..., 0) € K 9 , such that c T x = -A. With 
these specifications, the optimal LSD problem for rank-3 
states has been cast as a SDP. 

We proceed to show that this is a strictly feasible prob- 
lem. The state p has three positive eigenvalues and can 
be regarded as positive definite when considering only the 
subspace orthogonal to 7. We choose x T = (a, 0, . . . , 0) 
where < a/3 < smallest positive eigenvalue of p, so 
that 

1 

Pscp = a-13 > 0, 

1 , , 

p- Pscp = p-a-l 3 > 0. (24) 

The first and third blocks of F(x) are thus positive def- 
inite. For the second block, we need the fact that the 
eigenvalues of 7 Tl are given by \(1 ± p) and ±^q [2C| . 
Since < q < 1 and Q < p < 1, this means that 
Pscp — a |(l4 — 7 Tl ) is positive definite. Let us note 
that pjjp has zero eigenvalues only if q = 0, i.e., when 7 
is a product state. Thus, if we assume that 7 is not a 
product state, p^L > and we have strict primal feasi- 
bility. The case where 7 is a product state is treated in 
Sec. EE 



2. The dual problem 

The dual variable Z is now a 10 x 10 positive 
semidefinite matrix, subject to nine dual constraints. 
Strict dual feasibility is immediate as we can choose 
Z = diag(l3, 14,313), which can be easily checked to sat- 
isfy the nine dual constraints. 



3. Generalized Wellens-Kus equations 

Now, having established strict primal and dual feasibil- 
ity, we can invoke the complementary slackness condition 
([12]), or rather its rank-3 analog. The ith dual constraint 
now reads 

tr{F,Z} = itr^Zx + l 3 Z 2 Tl l 3 - Z 3 )} 

= a, i = l,...,9. (25) 

Any hermitian operator orthogonal to 7 can be written 
as i/rank3 = r i tr{TiH Yank3 } /ti{T 2 } . Let us repeat 

here that both Z 2 and Zj 1 have support in the total 
Hilbert space. To avoid inconsistency in the notation, let 
us define Zj^ , the restriction of the projection P 3 Z 2 1 P 3 

to its own support. We then arrive at Z 3 = Zi+Zj^ + 1L 3 . 
The third block equation in (I) of Eq. (fT"3")) then states 
that {Z\ + ZXl + t 3 )g pulc = 0, and since ^p Ure resides in 
the subspace that I3 projects onto, 

{Zi + -Z 2 |^)gp Urc = (?purei (26) 

and as before, this is supplemented by the eigenstate 
equation for Z 2 , 

(p-Q P urc) Tl Z 2 = 0. (27) 

Equations (f2"6"| and (f2"T)) are the generalization of the 
Wellens-Kus equations to the rank-3 case where the or- 
thogonal state is entangled. These are almost identical 
to the original equations, the subtle difference being that 
not only Z 2 , but also Z 2 ,} , the projection of its partial 
transpose onto the support of p, are now relevant. Sim- 
ilarly to the full-rank case, one can define Z\ + Z^ as 
the optimal entanglement witness for the state p. 

B. 7 is a product state 

1. The primal problem 

A little more care is needed if p is orthogonal to a 
pure product state 7 = |(l4 + (Ti)i(l4 — ti), the q = 
version of Eq. ([22]) . In this case, since p sep is separa- 
ble and orthogonal to 7, and 7 Tl must be orthog- 
onal too. The separability of p sep then requires: (i) the 
positivity of p^,p, and (ii) the orthogonality of p^_ and 
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7 Tl . Only two of the nine I\s do not obey Tj 1 ^/ 1 - 1 = 
These are T 8 = \{a 2 r 2 - a 3 T 3 ) and T 9 = ^(a 2 T 3 + <J 3 t 2 ) 
Furthermore, there exists a proportionality relation be- 
^ 7 Tl and r^ 7 

r 8 il 7 Xl = ^(r 8 il +ir^) = 

reads 



tween the products rJ 1 7 Tl and T^ 1 ^ 1 , inasmuch as 
'Ti-Ti _ i( F Ti -fiTj 1 ) = iT^ 1 7 Tl . Constraint (ii) then 



r sep / 



(x 8 - ixg)^ 1 



ir^) = o. 



(28) 



Since x% and xg are real, they must vanish and we have 
the parameterization p scp = ~ X)J=i x i^i- Consequently, 
the primal objective is now c^x with x G M 7 . The same 
choice of p scp = a^TL 3 shows that this modified primal 
problem is strictly feasible. 



aTg + for 9 appear as 



aT 8 + bT 9 = 



(31) 



;(|00> + |H)) 



where the two basis states are \<t> + ) = -^v 
and = -^(|01) + |10>). Equation ([29]) imposes 

that this matrix has an eigenvalue —1. Putting this re- 
quirement into its characteristic equation leads to the 
relation a + Ir = 1, and an angular parameterization, 
a = cos 0,6 = sin0, can be used. The corresponding 
eigenstate is nondegenerate, and can hence be identified 
with fjpm-o- Explicitly, we have 



2. The dual problem 

One can also verify, in the now familiar manner, that 
Z = diag(l3, iJ^Sla) is a strictly feasible point for the 
modified dual problem. 



\Q P u 



(32) 



which is maximally entangled. 

Now, Q scp has no components along r 8 and Tg, so one 
must have tr{E;(p — p pU rc)} = for i = 8,9. These turn 
out to provide a simple set of equations for the unknowns 
9 and S: 



3. Generalized Wellens-Kus equations 

The seven dual constraints lead to 
Z\ + Zj^ — Z3 + aT$ + bTg = — 13, where a and b 
are some real coefficients. We then arrive at another 
pair of generalized Wellens-Kus equations, 

(Zi + + ar 8 + bT 9 )g pul . c = ~g pluc , (29) 
(p-£ P urc) Tl Z 2 =0, (30) 

which are necessary and sufficient for optimality. Note 
that Z2 lies in the support of ij 1 while Zj 1 can have 
support in the total Hilbert space since Z 2 is not separa- 
ble. The term in parentheses in Eq. (f2"9"]) is the optimal 
entanglement witness for p. In contrast with the ear- 
lier cases, nonpositivity is provided by the combination 

zl\ + a r 8 + br g . 



4. Qsep has rank 3 

In the full-rank case, when the separable part is full- 
rank, the nonseparable part is maximally entangled. A 
similar property exists for rank-3 states. In three di- 
mensions, the analog of the full-rank case is a rank- 
3 state orthogonal to a pure product state. Note 
that the pure state has to be a product state to en- 
sure that all the relevant positive operators remain of 
rank 3 under partial transposition. If g sep has rank 
3, the optimal decomposition can be obtained analyti- 
cally. In this case, Z\ = Z 2 = and Eq. (|2"9")l reduces 
to (ar 8 + bTg)g pulc — —g pul - c . In the magic basis, which 
is a basis of Bell states, the nonzero matrix elements of 



tr{r 8 p} =(S-l)cos0, ( . 

tr{r 9 p} = (<S-l)sin0. 

Therefore, we obtain 5 = 1- ^/(tr{T 8 p}) 2 + {tr{T 9 p}) 2 . 
The solution to Eq. (|33[) then gives us g pmc and 
Qscp = P — £>purc in the optimal LSD of p. 

In general, one can assume that g sep has rank 3 and 
use the above result to determine the optimal g sep and 
£ purc . It is however necessary to check if the deduced 
g sep is indeed of rank 3 and separable. If the verification 
fails, g sep has rank 2 and one has to solve the generalized 
Wellens-Kus equations given in Eqs. - (I30|) . 



VI. CONCLUSION 

We have demonstrated that the problem of finding the 
optimal LSD of a 2-qubit state is a SDP. Indeed, the 
Peres-Horodecki criterion has permitted us to advanta- 
geously rephrase a separability constraint as a positivity 
constraint. We have shown that both the primal and the 
associated dual programs are strictly feasible, leading us 
to necessary and sufficient optimality conditions for LSD. 
In particular we have derived the original Wellens-Kus 
equations for full-rank states in a simple and natural way. 
Moreover we have generalized them to rank-3 states. We 
have also described the link between the dual SDP vari- 
ables and entanglement witnesses. Finally, many efficient 
algorithms for solving SDPs are available, allowing one 
to handle this problem numerically. Because the Peres- 
Horodecki criterion is also necessary and sufficient for 
composite systems of dimensions 2 x 3, it might be pos- 
sible to extend the SDP formulation to this case. 
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